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The Kondo-efFect is a many-body phenomenon arising due to conduction elec- 
trons scattering off a locahzed spin^. Coherent spin-flip scattering oflF such 
a quantum impurity correlates the conduction electrons and at low temper- 
ature this leads to a zero-bias conductance anomaly^*^. This has become a 
common signature in bias-spectroscopy of single-electron transistors, observed 
in GaAs quantum dot^a^a ^' — i ^ i ^' as well as in various single-molecule transis- 
tor a^°i^^i^^i^^i^^i^^ . While the zero-bias Kondo eflFect is well established it remains 
uncertain to what extent Kondo correlations persist in non-equilibrium situa- 
tions where inelastic processes induce decoherence. Here we report on a pro- 
nounced conductance peak observed at flnite bias-voltage in a carbon nanotube 
quantum dot in the spin singlet ground state. We explain this finite-bias con- 
ductance anomaly by a nonequilibrium Kondo-eff^ect involving excitations into 
a spin triplet state. Excellent agreement between calculated and measured non- 
linear conductance is obtained, thus strongly supporting the correlated nature 
of this nonequilibrium resonance. 

For quantum dots accommodating an odd number of electrons, a suppression of charge- 
fluctuations in the Coulomb blockade regime leads to a local spin-1/2 degree of freedom and, 
when temperature is lowered through a characteristic Kondo-temperature Tk, the Kondo- 
effect shows up as a zero-bias peak in the differential conductance. In a dot with an even 
number of electrons, the two electrons residing in the highest occupied level may either 
form a singlet or promote one electron to the next level to form a triplet, depending on the 
relative magnitude of the level splitting 6 and the ferromagnetic intradot exchange energy 
J. For J > 6, the triplet state prevails and gives rise to a zero-bias Kondo peak^, but 
when 6 > J the singlet state is lower in energy and no Kondo effect is expected in the linear 
conductance. Nevertheless, spin-flip tunneling becomes viable when the applied bias is large 
enough to induce transitions from singlet to triplet state. Such inter-lead exchange-tunneling 
may give rise to Kondo correlations and concomitant conductance peaks near V ~ ±6/e. 
However, since the tunneling involves excited states with a rather limited life-time, the 
question remains to what extend the coherence of such inelastic spin-flips and hence the 
Kondo-effect is maintained? 

In the context of a double-dot system, a qualitative description of a finite-bias Kondo- 
effect, leaving out decoherence and nonequilibrium effects, was given already in Reflld and 
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finite bias conductance peaks have already been observed in carbon nanotubesii^^^^ as well 
as in GaAs quantum dota^iSiiS. However, for lack of a quantitative theory for this nonequilib- 
rium resonance no characterization of the phenomenon has yet been possible. As pointed out 
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a bias induced population of the excited state (here the triplet), may change 
a simple finite-bias cotunneling step into a cusp in the nonlinear conductance. Therefore, 
in order to quantify the strength of correlations involved in such a finite-bias conductance 
anomaly, a proper nonequilibrium treatment will be necessary. As we demonstrate below, 
the qualitative signature of Kondo-correlations is a finite bias conductance peak which is 
markedly sharper than the magnitude of the threshold bias. 

We have examined a quantum dot based on a single- wall carbon nanotube (see Fig. la). 
Electron transport measurements of the two-terminal differential conductance were carried 
out in a cryostat with a base electron temperature of Tgi ~ 80 mK and a magnetic field 
perpendicular to the nanotube axis. The low temperature characteristics of the device 
are seen from the density plot in Fig. lb, showing dl/dV as a function of source-drain 
voltage V and gate voltage Vg. The dominant blue regions of low conductance are caused by 
Coulomb blockade (CB) while the sloping white and red lines are edges of the CB diamonds, 
where the blockade is overcome by the finite source-drain bias. Moreover, white and red 
horizontal ridges of high conductance around zero bias are seen. These ridges occur in an 
alternating manner, for every second electron added to the nanotube dot. These are the 
Kondo resonances induced by the finite electron spin S* = 1/2 existing for an odd number 
of electrons where an unpaired electron is localized on the tube. The zero bias resonances 
are absent for the other regions (with even A^) where the ground state spin is 5* = 0. Over 
most of the measured gate-voltage range (not all shown) the diamond-plot exhibits a clear 
four-electron periodicity, consistent with the consecutive filling of two non-degenerate sub- 
bands, corresponding to the two different sub-lattices of the rolled up graphene sheet, within 
each shellii. This shell-filling scheme is illustrated in Fig. Ic. 

We observe inelastic cotunneling features for all even A^ at eV ~ ±A for a filled shell 
and eV ~ ±5 for a half-filled shell. Reading off the addition, and excitation energies 
throughout the quartet near Vg = —4.90 V, shown in Fig. lb, we can estimate the relevant 
energy-scales within a constant interaction mode l^^i^^i^^ . We deduce a charging energy Ec ~ 
3.0 meV, a level-spacing A ^ 4.6 meV, a subband mismatch 5 ~ 1.5 meV, and rather weak 
intradot exchange, and intra-orbital Coulomb energies J, dU < 0.055. The fact that 6 > J is 
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consistent with a singlet ground-state for a half-filled shell, involving only the lower sub-band 
(orbital), together with a triplet at excitation energy 6 ~ J and another singlet at energy 
6. Notice that also the regions with odd and doublet ground-state exhibit an inelastic 
resonance at an energy close to 6. We ascribe these to (possibly Kondo-enhanced) transitions 
exciting the valence electron from orbital 1 to orbital 2, but we shall not investigate these 
resonances in detail. Focusing on the half-filled shell, i.e. N = 2, Fig. 2a shows the measured 
line-shape, dl/dV vs. V, at Vg = —4.90 V and for Tei ranging from 81 mK to 687 mK. 
The conductance is highly asymmetric in bias- voltage and exhibits a pronounced peak near 
V ^ 6/e which increases markedly when lowering the temperature (cf. inset Fig. 2b). 

For = 2, we model the nanotube quantum dot by a two-orbital Anderson impurity 
occupied by two electrons coupled to the leads via four different tunneling amplitudes tia, 
between orbitals z = 1,2 and leads a = L,R. In terms of the conduction electron den- 
sity of states, Up, which we assume to be equal for the two leads, the tunneling induces a 
level-broadening Fj = ^TTUpJ^a l^jaP- the Kondo-regime, Fj -C Ec, charge fluctuations 
are strongly suppressed and the Anderson model becomes equivalent to a Kondo-model in- 
corporating second order super-exchange (cotunneling) processes, like the one illustrated in 
Fig. 3a, in an effective exchange-tunneling interaction (see Supplementary Information, sec- 
tion 2). This effective model describes the possible transitions between the five lowest-lying 
two-electron states on the nanotube (see Fig. 3b). Since J <^ 6, we neglect J altogether and 
assume that the four excited states are all degenerate with excitation-energy 6. The next 
excited state is a singlet entirely within orbital 2 with excitation energy 26 (see Supplemen- 
tary Information, section 2). Although this state is merely a factor of 2 higher in energy 
than the five lowest lying states, we shall neglect this and all higher lying states altogether. 
Since there are no spin-flip transitions connecting the ground-state singlet to the excited 
singlet at energy 26, we expect this omission to have a minute influence on the conductance 
for V < 26. 

In order to deal with the logarithmic singularities arising in low-order perturbation the- 
ory within the Kondo model, it is convenient to employ the so-called poor man's scaling 
methodic. As we have demonstrated earlier^^, this method can be generalized to deal also 
with non-degenerate spin-states, like a spin-1/2 in a magnetic field, as well as with nonequi- 
librium systems where the local spin degrees of freedom are out of thermal equilibrium with 
the conduction electrons in the leads. The present problem is complicated by the presence 
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of the two finite energy scales, eV and S, causing the renormalization of the couplings to 
involve more than just scattering near the conduction electron Fermi-surfaces. To accommo- 
date for this fact, one must allow for different renormalization at different energy scales, and 
the poor man's scaling method should therefore be generalized to produce a renormalization 
group (RG) flow of frequency dependent coupling functions. The details of this method have 
been presented in Refs. I25i26l and the specific RG-equations for this problem are presented 
in section 3 of the Supplementary Information. 

These coupled nonhnear RG-equations do not lend themselves to analytical solution but 
may be solved numerically with relatively little effort. From this solution, we obtain renor- 
malized coupling-functions with logarithmic peaks at certain frequencies determined by eV 
and 6 (see e.g. Fig. 3c). Since resonant spin- flips take place only via the excited states 
at energy 6, Korringa-like spin-relaxation, via excitation of particle-hole pairs in the leads, 
partly degrades the coherence required for the Kondo-effecli^i. Even at zero temperature, 
the couplings therefore remain weak and increase roug hly as l/\og{VWTT^/TK), with F 
being the effective (\^-dependent) spin-relaxation rate. Using the parameters tai from the 
fit in Fig. 4, we calculate the Kondo-temperature and the spin-relaxation rate and find that 
Tk = 0.4 mK and F ^ 350 mK (see Supplementary Information, sections 3-5). This in 
turn implies a small parameter 1/ ln[5/Tft:] ^0.1, which justifies our perturbative approach. 
Notice that the smallness of corresponds roughly to a mere 50% reduction of the total 
hybridization to orbital 1 as an electron is added to the nanotube. This tiny Tk does not 
show up directly in the V^-dependence of the nonlinear conductance, which is characterized 
instead by the excitation energy 6 and the spin-relaxation rate F. 

Having determined the renormalized coupling functions, the electrical current is finally 
calculated from the golden rule expression: 
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(V^-V), (1) 



where are the nonequilibrium occupation numbers of the five impurity states and 
dJia'-Lai^) renormalized exchange-tunneling amplitude for transferring a conduc- 

tion electron from left to right lead and changing its spin from cr to cr', while switching 
the impurity-state from 7 to 7'. Figure 4 shows a comparison of this calculation to the 



data, obtained by fitting the four tunneling amplitudes Ua determining the unrenormalized 
exchange couplings. Notice that in practice the measured conductance at V — and at 
\V\ 3> S, together with the asymmetry between positive and negative bias, places strong 
constraints on these four amplitudes. The resulting fit is highly satisfactory, with only slight 
deviations at highest voltages. 

To quantify the importance of Kondo-correlations, we compare also to plain nonequi- 
librium cotunneling. This unrenormalized second order tunneling mechanism clearly 
underestimates the conductance peak at positive bias, and most severely so at the lowest 
temperature. Qualitatively, simple nonequilibrium cotunneling gives rise to a cusp in the 
conductance which is roughly as wide as the magnitude of the threshold-bias itself. As 
evident from the data and analysis presented here, nonequilibrium Kondo-corrclations 
instead produce a conductance anomaly which is sharper than the magnitude of the 
threshold,- that is a proper finite-bias peak. 
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Figure 1 Experimental setup and shell-filling scheme for a single-wall carbon 
nanotube. a, Schematic of the nanotube device, comprising a single-wall carbon nanotube 
grown by chemical vapor deposition on a Si02 substrate and contacted by Cr/Au source and 
drain electrodes, spaced by 250 nm. Highly doped silicon below the Si02 cap layer acted as a 
back gate electrode. Room temperature measurements of conductance as a function of back- 
gate voltage, Vg, indicate that the conducting nanotube is metallic with a small gap outside 
the region considered here (cf. Supplementary Information, section 1). b. Density-plot of 
dl/dV as a function of V and Vg at Tgi = 81 mK. c. Diagram illustrating the corresponding 
level-filling-scheme with defined as the total number of electrons modulo 4. 
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Figure 2 Temperature and magnetic field dependence of the finite-bias reso- 
nance, a, dl/dV as a function of V, at Vg — —4.90 V (dashed cross section in Fig. lb) , taken 
at Tel =81 (thick), 199, 335, 388, 488, 614, 687 (dashed) mK. b. Temperature dependence 
of the finite-bias conductance peak in Fig. 2a (open circles) and the neighboring zero-bias 
Kondo resonance (solid circles) a.tVg = —4.96 V (see Fig. la). The open circles are consistent 
with the saturation of a log-enhanced finite-bias Kondo-peak as T becomes smaller than the 
spin-relaxation rate F ~ 350 mK, determined from the parameters in Fig. 4. The solid 
circles follow the NRG-interpolation-formula a(l + (2V0-22 _ i)(r/r^)2)-o.22 ^^^^ ^ ^ g.ll 
and Tk = 1.0 K (dashed fine), c, dl/dV vs. V and magnetic field B, at Vg = -4.90 V. d, 
dl/dV vs. V ai B — 6 T, corresponding to vertical cross section (dashed) in Fig. 2c. The 
d'^I/dV'^ trace (thin) underhnes the presence of three distinct peaks in dl/dV. The black 
bars each correspond to AV = g/iBB/e with g = 2.0 for nanotubes. 
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Figure 3 Schematic of the inelastic spin-exchange underlying the Kondo-effect. 

a, Illustration of the cotunneling mechanism giving rise to the effective exchange-interaction 
in the Kondo model. Note that the virtual intermediate state (red spins) in panel 2 is 
suppressed by a large energy-denominator of the order of the electrostatic charging energy 
of the nanotube. b, Schematic of the five different low-energy states retained in the effective 
Kondo-model: \s) is the singlet ground-state, |m) (m=l,2,3) are the triplet components with 
excitation energy 6 — J ^ 6 and \s') is the corresponding singlet with excitation energy 5. 
Applying a magnetic field causes the triplet to split up as seen in Fig. 2c, d. c, Renormalized 
dimensionless exchange-coupling vs. incoming conduction electron energy u, measuring the 
amplitude for exchange-tunneling from right to left lead while de-exciting an electron on the 
nanotube from orbital 2 to 1. Different colours encode the different stages of the RG-flow 
as the bandwidth is reduced from Dq to zero. Parameters as in Fig. 4. 
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Figure 4 Fitting the nonlinear conductance by perturbative RG calculation, a, 

Experimental data (red dots) at the lowest temperature, fitted by the perturbative RG 
calculation (white dots) with {^li, tL2, tm, ti?2} = {0.024, 0.033, 0.111, 0.083} v^^c/z/f and 
T = Tel = 81 mK. We use J = and an arbitrary bandwidth, Dq = 1 eV, larger than all 
other energy-scales. From the bias at which dl/dV vs. V has maximum slope, we read off 
S = 1.515 meV. The blue curve represents simple cotunneling, i.e. second order perturbation 
theory(PT) with tunneling-amplitudes tia fitted to the data at T = 687 mK (see panel b). 
b. Comparing to the data (red and white dots) at the highest temperatures, using same 
parameters as in panel a, but with T = = 687 mK. The blue curve represents the best 
fit of second order perturbation theory with T = 687 mK. 
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Supplementary Discussion 

1. Nanotube transport experiments 

Our nanotubes were grown by chemical vapor deposition (CVD) on a doped silicon sub- 
strate with a 400 nm oxide cap layer. Ferric iron nitrate nanoparticles deposited from a 
solution in isopropyl alcohol acted as catalyst for the CVD process, which was carried out in 
a tube furnace by flowing methane and hydrogen over the sample at 900 °C— . This process 
yielded mostly individual single-wall nanotubes as determined by atomic force microscopy. 
The nanotubes were contacted by thermally evaporated metal electrodes (35 nm Au on 
4 nm Cr), spaced by 250 nm and patterned by electron beam lithography. The two-terminal 
conductance was measured using standard lock-in techniques with ~ 5 ac excitation 
and voltage bias V applied to the source with the drain grounded through a low-impedance 
current amplifier. The effective electron temperature can differ from the measured mixing 
chamber temperature in the dilution refrigerator. A base electron temperature of T^i ~ 
80 mK was estimated from the temperature dependence of device characteristics. 

In the range of gate voltage considered in this paper, the room temperature conductance 
of the device is around 1.8 /h and independent of back-gate voltage Vg. At lower gate 
voltages, a weak Vg dependence at room temperature is seen and can be attributed to a 
small band gap in the metallic nanotube induced by perturbations.— The coupling decreases 
with Vg and at low T the Coulomb blockade peaks disappear at ~ —2.5 V, i.e. around 
AA^ = 40 electrons away from the present Vg range.-^ 
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2. Low energy two-particle states and effective Kondo-Hamiltonian 

The relevant single-particle states of the nanotube carry spin a =t, J, and orbital index 
i = 1, 2, and the lowest lying two-particle states are denoted as follows: 



\s) =(|i,T)-|T,i))®|i,i)/V2 (2) 

-1) =|i,i)®(|l,2)-|2,l))/v^ (3) 

|0) =(|i,t) + |t,i))®(|l,2)-|2,l))/2 (4) 

|1) =1 T,T)®(|1,2)-|2,1))/V2 (5) 

\s) =(U,T)-|T,i))®(|l,2) + |2,l))/2 (6) 

\s") =(|i,T)-|t,i))®|2,2)/V2 (7) 



These states have energies E^^ i?_i,o,i = -Es + 5 — J, E^i = Es + 6 and Eg" = Es + 26, and 
since J < S, the singlet | s) is the ground-state. We include only the five lowest lying states 
and neglect the highest lying singlet, | s"), altogether. Within the Hilbert-space of these five 
states, the effective Kondo-Hamiltonian takes the form: 

k,(T 
1=1,2 
cc=L,R 

where 



with 



k,k',cr,cr' 
a,a'=L,R 



m=— 1,0,1, s' 



9a'a = '^^Ftia't*jEc, 



ail = {ali-r = 2^vpt,^diJEc, 

„12 ^ / 21 ^* = ^12 
fa' a xfaa' I Ha'ai 



0. 



(8) 



(9) 



The vector of Pauli-matrices is denoted by Tij and all terms in the interaction part have the 
form of spin, and orbital exchange, except for the two terms in the last line proportional 
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to SijSa'a which are pure potential scattering terms. Throughout, we use the convention 
A± = Ax ^ iAy, for any vector-operator A. The Hamiltonian is expressed in terms of the 
two-particle vector-operators 

5+ = (5-)t = v^(|l)(0| + |0)(-l|), 

= |i)(i|-|-i)(-i|, 

P^2 = {P^i^ = V2\l){s\, 

Pi2 = iP2iy = -\o){si 

T+ = (T-)t = V2(-|l)(s'| + |s')(-l|), 
= |0)(.'| + |s')(0|, 

together with the scalar-operators 

\s){s\ = liP'-ls% 

\s){s'\ = (|.')(s|)t = -ip.f, 

E \m){m\ = ^iT' + S% (10) 

m=-l,0,l,s' 

where P = Pu + Ai- Defining also the operator M = [T,P]/3i = i{\s){s'\ — \s'){s\), we 
note that the ten operators {S, T, P, M} satisfy the commutation relations: 

[Si, Sj] = iSijkSk, [Si, Tj] = lEijkTk, [M, Tj] = iPi, 

[Ti,Tj] = ieijkSk, [Si,Pj]=ieijkPk, [M,P,] = -iTi, (11) 

[P„ P,] = te,,kSk, [r„ P,] = i6,,M, [M, Si] = 0. 

which identifies them as generators of the Lie-algebra S0{5). Notice that T'^+P'^+S'^+M'^ = 
4 is a Casimir-operator, i.e. a constant of motion. The Hamiltonian (jH} is manifestly 
invariant under spatial rotations, but it will posses this abstract S'0(5)-symmetry only 
when 6 = J = and tia = t2a, in which case one would expect to observe a conventional 
zero-bias Kondo-peak, characterized by a Kondo-temperature which is enhanced compared 
to the usual SU{2) Kondo-effect. Notice that an S'0(5)-Kondo-effect has been discussed 
earlier in the context of a triple-quantum-dot systent^. 



— # — * 

As pointed out for a double-dot system studied in Ref. |32, 5" and P generate SO (4), 



and it is the addition of the excited singlet \s') which adds four new generators in the 
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present problem. In a double-dot system at even filling, the inter-dot tunneling breaks 
the degeneracy between singlet and triplet states and an inelastic cotunneling channel 
is generally available. A finite-bias Kondo-resonance in such a double-dot system was 
suggested in Ref. Q and the current due to nonequilibrium cotunneling was calculated in 
Ref. |2l|, but the combined nonequilibrium Kondo-effect has not yet been examined. By 
simply leaving out the excited singlet |s'), our present calculation could readily be applied 
to this problem as well. 



3. Perturbative renormalization group equations 

The (one-loop) perturbative renormalization group (RG) equations satisfied by the fre- 
quency dependent couplings are established in much the same way as explained earlier in 



Ref. 



2d Here, they take the following form: 



d\nD 



dlnD 



dlnD 
dhiD 



dlnD 



Y.{9%Ac^"Vl2)g%^{a"Vm 



uj-a"V/2 



1 



9X"{<^"V/2) + - (pL^a"(«'V/2) - q:,Ac^"V/2)) 



+ 9l'M"V/2) - ]- U"M"V/2) - C,M"v/2) 



-\ff:}^9tA(^''V/2)C,M''V/2 - 46)e^^^^S-a"V/2: 
a" 

-I E { 59A"ic^"V/2)g%M"V/2 + 4S)0u.-r^5-a"V/2 
a" 

+39tA(^''V/2)9A{a''V/2)e^.^.v/2 
+tIpI.Ao^"V/2)p%M"V/2 + Tl5)e^_,.^s-."v/2 
- rM."(«'V/2)pJ,„(a'V/2)e^.„„v/2} , 

\rl J2 { 39tAo^"V/2)C,M"V/2 - t^5)^.+^5_„.w2 

a" 

+ piAo^"V/2)j}A{o^"V/2 - rl5)e^^,^^s-."v/2] , 



(12) 



(13) 
(14) 



(15) 



(16) 



with the shorthand notation 9^ = 9{D — |x|), 1 = 2 and 2 = 1. For a given set of initial 
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values (at scale D = Dq), parametrized by the bare tunneling amplitudes tia according to 
(ini), these equations are readily solved numerically for arbitrary uj and D. Taking the limit 
of 0, we obtain the renormalized coupling functions used in the Golden Rule expression 
for the current. 

As usual, when deriving the Kondo-model from an Anderson model, it is convenient to 



introduce two angles, (cos 0i, sin 0j) = {tn, Ur) / \/t1i^ + t^/j, and parameterize the exchange 
couplings as 

{da'a^ a > Ca) = idij , Pij , Qij} i . ^]{ COS (pj siu C^j ) (17) 



with initial conditions g,j = 2(z/i./Ec)y'(l + ri)(t-i + t-^)(tjx + ^ji?)) Pa = ^ii9iu P12 = 
P21 = 912 and Qij = 0. The L/R matrix-structure of the couplings is now exterior, in the 
sense that the RG-equations are identical for every (a', a) component. 

For D ^ 6, V,T,u, the RG-equations simplify to describe the flow of coupling constants 
{g,p,q)ij with no dependence on frequency, uj, nor lead-index, a = L,R. The difference 
in tunneling-strengths to the two orbitals still plays an important role, and enters via the 
initial conditions. From the RG-equations it is clear that all the different couplings diverge 
at the same energy-scale and it is this scale which we henceforth refer to as the Kondo- 
temperature, Tk (see Fig. E)). The numerical solution of these high-D RG-equations allows 
us to study the dependence of Tk on the three independent initial values gu, g22 and Dq. 
This reveals the following nearly perfect interpolation-formula for T^: 

where 1.6 is a good approximation to a universal number describing this particular Kondo- 
effect. For gu = g22 this reduces to Tg_~ Dq exp{—0.36 / gu) , as found earlier in a study of 
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vertical quantum dots^'^(cf. also Refs. |3j|35|). 

Notice that in the standard Kondo-effect involving a zero-bias conductance peak, the 
Kondo-temperature can be estimated directly from the width of the conductance peak. In 
the present problem, however, Tk is of little physical significance insofar as the non-linear 
conductance is characterized mainly by the spin-relaxation rate F (see sec. 5) and the 
subband-mismatch 6. Nevertheless, Tk still encodes a scaling-property of the RG-solutions 
and shows how our rather arbitrary choice of bare bandwidth Dq = 1 eV is linked to 
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FIG. 5: Inverse couplings vs. bandwidth D. The inverse couphngs all vanish at the 
same energy-scale defining Tr-. The couplings shown here are based on the same bare pa- 
rameters as were used for the low temperature fit in the main paper: {t^i, tL2, ii?2} = 
{0.024, 0.033, 0.111, 0.083} y^c/i^F, J = and = 1 eV, which implies that Tk « 0.4 mK. 

the values of the bare couplings. In spite of the fact that the Kondo-temperature in this 
problem is enhanced over that of the standard SU(2) symmetric Kondo-model, our fitting 
parameters imply that Tk ~ 0.4 mK, which is much smaller than the value of 1.0 K found 
for the neighboring Coulomb-blockade valley having an odd number of electrons on the tube 
(cf. Fig. 2b in main text). Estimating the effective coupling in the neighboring valley from 
Eq. p8|) with (722 = 0, we find that this big difference in Kondo-temperature corresponds to 
a mere reduction of the total hybridization to orbital 1 (a/^Ii + t^n) by roughly 50% when 
changing from N=l to N=2. 

4. Nonequilibrium distribution functions 

Using the renormalized coupling-functions, we may calculate the transition rates between 
the various two-particle states using the Golden Rule expression 

rj 1*00 

W,,,{V,5,T) = ^ / du;Y,\9t:%.,M?fi^-^^c.){l-f{uJ + e,-e,,-^^o.,)), (19) 

J -oo I . 

a,a'=L,R 

where / denotes the Fermi-function. The nonequilibrium distribution functions for the two- 
particle states, n-y, are then found by solving the steady-state quantum Boltzmann equation 
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y^^W^yny = y^^Wy^n^, (20) 

7' 7' 

together with the constraint ^y=s -i o i s' ''^i ~ ensuring that the two electrons in the 
half-filled shell on the nanotube occupy exactly one of the five lowest two-particle states. 
These volt age- dependent distribution functions are subsequently plugged into the Golden 
Rule expression for the current. 



5. Spin-relaxation 

As demonstrated in Ref. |27|, a finite bias gives rise to Korringa-like spin-relaxation via 
inter-lead particle- hole excitations and for a single spin-1/2 the logarithmic singularities were 
found to be cut off by the (V-dependent) spin-relaxation rates 1/Ti 2. To be more precise, 
it is the broadening of the transverse dynamical susceptibility, I/T2, which cuts off the log- 
renormalization of spin-ffip exchange couplings, whereas the renormalization of non-spin- 
ffip exchange couplings are contained by the broadening of the longitudinal susceptibility, 
1/Ti. Diagrammatically, these rates arise from a combination of both self-energy, and vertex 
corrections to the spin-susceptibility bubble and a simplification which simply omits the 
vertex corrections will lead to serious mistakes. 

In the present problem, the inter-orbital transition-rates (from excited to ground-state) 
are strongly enhanced by the large available phase-space. Ws,m = —1,0,1) and Ws^s' 
turn out to be larger than all other transition-rates by roughly an order of magnitude and 
the self-energy broadenings or line-widths of the excited states are therefore much larger 
than that of the ground-state. At least one of these large self-energy broadenings will 
contribute to the cut-off in the log-renormalization of all exchange couplings except one. 
Pa couples the ground-state to itself and the relevant self-energy broadening is very weak. 
For this particular coupling, vertex corrections will now play the dominant role and they 
supply additional terms in the total broadening involving the large Ws^m and Wg^s'- AH 
log- singularities are therefore cut off by rates of the order of these dominant line-widths of 
the excited states. 

From the expression for the current (cf. main text Eq.(l)), it is clear that the conduc- 
tance peak at ~ (5 is mainly determined by the inter-orbital exchange-couplings (7*' and 
p". Furthermore, these couplings turn out to be larger than the others and altogether the 
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main influence of spin-relaxation on the physical current is therefore via these inter-orbital 
couplings. In other words, broadening all log-renormalization by a single effective spin- 
relaxation rate estimated from an inter-orbital susceptibility is expected to produce only 
minute errors. 

In terms of a generalized ss'-susceptibility, we determine the physical spin-relaxation-rate 
Ts^s' as the sum of self-energy, and vertex corrections: 

r..' = rr' + liYl ^7',^ + E ^7',^'), (21) 

where F^'**' is the contribution from vertex corrections. Just as the vertex-corrections to 
the transverse spin-susceptibility for a single spin-1/2 include only non-spin-flip processes^I, 
the correction F^'*' picks up only intra-orbital transition-rates and is therefore negligible 
compared to the dominant term Ws,s' coming from the self-energy broadening of the excited 
state. Tsm = —1, 0, 1) is determined in a similar way and since these two different inter- 
orbital relaxation rates turn out to be very close in magnitude as well as in V-dependence 
we define a single effective spin-relaxation rate as their average: 

F = i(F,,, + F,o). (22) 
The spin-relaxation mechanism is incorporated in the RG-equations by replacing 6^ by 



9{D — I v^eM-T^+T^I) in Eqs. p2H16|) . Furthermore, all Fermi-functions occurring in 
the transition-rates and the current are effectively smeared by F, by replacing the energy- 
conserving (5-functions appearing in the Golden rule expressions with Lorentzians of width 
F (i.e. spin-susceptibilitites Im[x^]). For example, we evaluate the transition rate Wst as 
follows: 



Q /•OO /»00 

Wst{V,5,T) =-Y, / du^ rf^|^7S:>,.(^)r/(^-/^«)(l-/(c^ + ^-/Xa'))Im[xfT(^)] 




1 



F/tt 



{e - 5)2 + F2 

, . , J J -k fc -1- i 



h 

a,a'=L,R 
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where N is the Bose-function and A = \^P~+~J^~+V^~+T^ is an ultra-violet cut-off on 
the spin-susceptibility ensuring convergence of the integral over e. Notice that, in order to 
speed up the numerical evaluation, we have replaced the square of the coupling-function 
by its average over the window set by the Fermi-functions. The error introduced by this 
approximation is estimated to be subleading in the small parameter 1/ log(5/T/^). 

We find that r{V — 5) ~ 350 mK (varying between 345 mK and 475 mK when changing 
V over the measured range) and the data clearly sample the full crossover from low to high 
temperatures with Tjf^^«* 81 mK < T < 687 mK ^ T^^9hest_ 
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